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Abstract. 

The second Z 3 parafermionic conformal theories are associated with the coset con- 
struction su (2h ><su{2)4, ^ s iid_ on _ so iid integrable lattice models obtained by fusion of the 
model based on level- 1 representation of the affine algebra have a critical point 
described by these conformal theories. Explicit values for the Boltzmann weights are de- 
rived for these models, and it is shown that the Boltzmann weights can be made positive 
for a particular value of the spectral parameter, opening a way to eventual numerical 
simulations of these conformal field theories. Away from criticality, these lattice models 
describe an integrable, massive perturbation of the parafermionic conformal theory by 
the relevant field ,^Di 3 . 



Over the past twenty years Conformal Field Theories (CFT) have proven to be a very 
powerful tool to describe two dimensional statistical systems at criticality. Recently a 
new way to study critical systems has been developed : Schramm Loewner Evolution 
(SLE) [H[2]- In this approach to conformally invariant systems, clusters are described 
through the conformal properties of their domain walls. While the connection between 
minimal models (c < 1) and SLE is relatively well understood, it is still an open problem 
for CFTs possessing extra symmetries. Statistical systems with a discrete symmetry 
Z*n are believed to be described by parafermionic theories. The first of such theories, 
constructed by Fateev and Zamolodchikov [3J, have been studied in the context of SLE 
[SI [6] , and using particular lattice realizations of the first parafermions, numerical studies 
of fractal interfaces have been performed [3 [9] . 

Conformal theories associated with the coset construction S0 ^1^ T ^ ^ 2 (N > 5) are 
referred to as (second) parafermionic CFTs Z N (k) [HI [El ESJ E]. Having a lattice 
version of these parafermions is very appealing for several reasons. It would allow a 
better understanding of this specific Zjy symmetry, how it can be realized, and how it's 
broken when the phase transition occurs. And it gives access to the related SLEs, at 
least numerically through the study of interfaces at criticality. 

Using the coset description of the parafermions, techniques to build integrable lattice 
models are available : the fusion procedure [TS1 IS] applied to particular existing lattice 
models [17]. In general the models we end up with have fluctuations variables both on 
sites and edges of a square lattice, and local interactions are around faces. When N = 3 
this general picture is somewhat simpler, since there is no edge variable in that case. The 
corresponding restricted solid-on-solid (RSOS) models have been studied in the context 
of solvable models obtained by fusion of the 8 vertex model [THJ [THl ED] • In particular 
their critical behavior has been identified by evaluation of the local state probability 
(LSP), and it is known that they are described by the N = 3 second parafermions. 

(2) 

The purpose of this article is to study an explicit realization of these second 
parafermions [TU] . We adopt a somewhat different point of view, having in mind a 
generalization toward the case N > 3, but up to some conventions the lattice model we 
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get is strictly equivalent to those described above. The focus of the paper is on practical 
aspects of this model, in particular it is shown that all Boltzmann weights can be made 
positive for a particular value of the spectral parameter, opening a way to numerical 
simulations. 



1 Description of the model 

Starting from the coset description of the parafermions : Z^(fc) = > the 

standard way to build an integrable lattice model is by fusion of the model realizing 
S ° 50(JV)t+[ jV ' >1 " These JMO models [T7j are related to the vector representation of SO(N). 

For N = 3 it is known that this is equivalent to the model corresponding to the 
symmetric tensor representation of degree 2 of SU(2), which in turns is obtained by 
fusion of the model based on the vector representation of SU{2). This last model is 
nothing but the celebrated ABF model |22j . 



(2) 

There are two equivalent ways to obtain a solvable lattice models realizing Z 3 . by 
a 4- fusion of the vector SU(2) model, or by a 2-fusion of the vector 5*0(3) model. In 
terms of CFT, this is related to the following coset equivalence : 

(2) _ SU(2) k x SU(2) 4 _ 50(3) fc/2 x 50(3) 2 
3 { ) ~ SU(2) k+4 ~ SO(3) k/2+2 

As a coset based on SU(2), integrable lattice models obtained by 4 fusions of the 
ABF model are readily available [18j [20] . 

Equivalently, this can be seen as a 2 fusion of the model based on vector representation 
of 5*0(3), and this is the point of view adopted here, since it will be closer to the general 
case (N > 3). This will lead to some different conventions from the usual SU(2) lattice 
models, in particular we will have to allow half integer values for the heights and for the 
"level" k/2, but it should be stressed that both models are exactly the same. 

Consider a two dimensional square lattice, with heights a, b, c, d, etc at each site 
taking the values 1/2, 1, L — 1/2. These local states are subject to nearest neighbour 
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Figure 1: Nearest neighbour constraints for L = 17/2. Only the integer heights are 
shown. Arrows indicate admissible pairs. 
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Figure 2: Four sites surrounding a face. The sublattices A and B are shown. 

constraints : if a and b live on neighbouring sites, the pair (a, b) must be admissible, 
i.e. a — b = ±2, ±1, and 2 < a + b < 2L — 2 (see Fig. [I] for an exemple). For this 
model to have any internal freedom, the parameter L must be large enough : L > 7/2. 
In fact L will be related to k by L = k/2 + 3, so that the tricritical three state potts 
model correspond precisely to L = 7/2. It turns out that due to the nearest neighbour 
constraints, for a given configuration the heights will be either all integer, or all half 
integer. We end up with two models for a given value of L, depending on whether we 
consider integer or half integer values for the heights. 

The restricted solid-on-solid model is a face model, in the sense that the heights 

variables interact around faces of the square lattice. The statistical weight assigned to 

an elementary face of the lattice is zero unless all four pairs of adjacent heights on the 

( a b\ 

edges are admissible. We will denote these weights as W \ where a, 6, c, d are 

\d c J 

the four surrounding sites, ordered clockwise from the upper left, as in Fig. [21 

Such a model is insensitive to the following gauge transformation of the weights, since 
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the additionnal edge interactions will cancel out on two neighbouring faces 



W 




f(a,b)g(a,d) 



f(d, c)g(b, c) 



W 




(1.2) 



To be more precise, the partition function on a torus (i.e. with periodic boundary 
conditions) is unchanged. If fixed boundary conditions are applied on the boundary, the 
partition function is simply multiplied by an irrelevant factor. This is typically the case 
when calculating the local state probability. If one is interested in more general boundary 
conditions, such as open boundary conditions, then these gauge transformations have a 
non trivial effect, localized at the boundary of the system. In the bulk, the model is 
insensitive to these gauge transformations. 

Such transformations will be required to ensure the positivity of the Boltzmann 
weights (cf section [2]). 

It should be noted that these gauge transformations preserve integrability : if the 
weights obey the star-triangle equation, they will still satisfy it after transformation. 

Explicit values for the Boltzmann weights are obtained by fusion procedure of the 
original SO(3) model. These are given in the appendix lAl 



exactly solvable models, they are usually called spectral parameter, and nome (respec- 
tively). Whereas the spectral parameter u plays a crucial role in the integrability of the 
model, it can be fixed once and for all for numerical simulations. The nome < p < 1 
plays the role of temperature. The system is critical for p = 0, and at zero temperature 
for p = 1. We will consider here only the so-called regime III, where — 1/2 < u < and 
< p < 1. 

By construction, the fusion procedure preserves many interesting properties from the 
original model, the most important of them being integrabililty : the weights of the fused 
model obey the star-triangle relation (STR) Fig. [3J The STR ensures that, in the case 
of periodic boundary conditions, row-to-row transfer matrices for different value of the 




depend on two parameters u,p. In the context of 
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Figure 3: The star-triangle equation. The no me p does not a play a role in the STR and 
has been omitted : it is the same for all faces. 

spectral parameter u commute, the nome p being fixed. 
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(1.3) 



The Boltzmann weights obtained by fusion enjoy the following properties, herited 
from the original model : 

• Reflection symmetry : 
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Rotational symmetry 
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• Z 2 symmetry : 
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(1.6) 



where the following notations are used : A = —1/2, G a = [a], and [x] denotes an 
elliptic theta function (cf Appendix 1X1) . 

From the rotational symmetry it appears that the spectral parameter u is related to 
the spatial anisotropy of the interaction. The model is istropic for u = A/2 = —1/4, 

1 /2 

since the additionnal terms ( < ^ b< ^ d j in equation (jl.5p can be absorbed in a gauge 
transformation (cf section [2]). 



2 Some useful gauges transformations 

From the values in the appendix, it appears that some Boltzmann weights are not pos- 
itive at u — —A/2. But in fact these signs can all be absorbed in the following gauge 
transformation : 



w\ ab \ %^4w l a " ] (2.1; 

d c J e(d,c)e(b,c) y d c 

where e(a, b) = e(6, a) is just a sign, given by (we consider here the case when the 
heights are integers) : 



e( a ,a) = (-1)" 
e(a,a+l) = (-i)^ 1 )/ 2 

e(a,a + 2) = 1 (2.2) 

This works with half integer just by shifting the heights by 1/2, defining e(a,b) = 
e(a — 1/2, b — 1/2) when a, b are half integers. 
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The only effect of this gauge transformation is to change the sign of the following 
weights, and their symmetric obtained by reflexion and rotation : 



W 



a a 
a± 1 a±2 
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(2.3) 



In particular, this gauge preserve all the symmetries (reflexion, rotation, Z 2 ) of the 
Boltzmann weights, and ensures their positivity for u = —1/4. The model is then both 
isotropic and positive for this specific value of the spectral parameter. Finally, to make the 
model explicitly isotropic, the following gauge transformation, with f(a, b) = {G a / 'Gb) 1 ^ 4 \ 
g(a,b) = 1, can be performed : 
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(2.4) 



This preserve Z 2 and reflexion symmetries, and positivity, since G a = [a] > for all 
heights. And it makes the rotational symmetry more explicit : 
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(2.5) 



Once the spectral parameter is fixed at the isotropic value u = —1/4, the only re- 
maining parameters of the model are the nome p, < p < 1, and the boundary on the 
heights L. If one is interested in the critical behavior of the model, the nome p must 
be put to 0, in which case the Boltzmann simplify considerably : all theta functions [x] 
merely become sin(7nr/L). To sum up, the model is critical and isotropic for u = —1/4 
and p = 0. The only parameter is then L = k/2 + 3, and correspond to the choice of the 
CFT describing the continuum limit of this model : Z^\k), k > 1. 
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3 Low temperature limit 



The low temperature limit correspond to p — > 1. The following reparametrization of the 
elliptic theta function is very useful to study this limit : in regime III, the weights can 
be reexpressed using the parameter x rather than p, where : 



p = e ~ e x = e~ 4n/Lt 



p — > 1 x — > 



To do so we use the following modular transformation : 



(a-L/2) 2 r 

[a] = t{x)x~^ E(x a ,x L ) 

oo 

E(z, x) = Jl (1 - zx 11 - 1 ) (1 - z- l x n ) (1 - x n ) 

n=l 



where r(x) is a certain function of x here, irrelevant here because it cancels out when 
evaluating ratios of theta functions. The zero temperature limit is now x — > 0. This 
is very convenient to study the low temperature properties of the model, because the 
product expansion for E(z,x) is rapidly convergent as x — > 0. 

In particular the limit x — > 0, w = x u fixed, is very interesting since the weights 
become diagonal : 



lim W 

x^0,w=cst 



I a 
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u,p^j 
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u,p) = F(u)5 b4 w- H{aM+fa+fc - fb - fd 



where 



F(u) 

fa 



X L 



\a-L/2\' 
2L 



In this expression F(u) is just a global multiplicative factor, and w^ a+ ^ c ' h * A can be 
absorbed in a gauge transformation. What remains is : 



In this limit, only configurations invariant under translations along the southwest to 
northeast diagonal contribute to the partition function. Thus what is left is an effective 
one-dimensional problem. Using the corner transfer matrix trick, one can use this one- 
dimensional limit to calculate exactly one point functions (cf appendix [C]). Let us proceed 
to examine the ground states of this model. 

Ground states are configurations that maximize the Boltzmann weights (in absolute 
value). In regime III, which correspond to the region — 1/2 < w < 0, < p < 1, the 
ground states do not depend on the value of u and p. The one-dimensional limit is very 
practical to extract information about them. One can readily see that all ground states 
must be invariant under translations along the southwest to northeast diagonal. Then 
to maximize the term must take (a, b, c) such as H(a, b, c) is minimal : this 

is acheived for a = c. Ground states are thus invariant under translations along the 
southeast to northwest diagonal as well. 

This is consistent with the fact that for the regime III, ground states must be isotropic, 
since one could have been interested in the limit where x — > 0, while w = x x ~ u is fixed. 
In that limit the weights would behave as (up to some gauge transformation): 



Every configuration invariant under translations along diagonals is a ground state. 
Decomposing the square lattice in sublattices A and B (see Fig. [1]), a ground state will 
be described by a pair of admissible heights (b, b + r/), such that heights assume the value 
b on the sublattice A, and c = b + r\ on B. Thus there are five types of ground states : 




(3.1) 
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Figure 4: Ground state (6, c). Heights living on the sublattice A assume the value b, 
while those on the sublattice B take the value c = b + 7] 

• (6,6), for 3/2 < b < L - 3/2 

• (6, 6 + 1) and (6 + 1, 6) for 1 < 6 < L - 2 

• (6,6 + 2) and (6 + 2,6) for 1/2 < 6< L - 5/2 

The boundaries on 6 come from the admissibility rules for neighbouring sites. A 
compact way of writing these constraints is: 3/2<6<L — 3/2. where 6 is just the 
average height value in the ground state (6, c) : 6 = 2±£. It should be stressed that (6, c) 
and (c, 6) are different ground states (for c 7^ 6) . 

It is interesting to note that there is no obvious realization of the Z3 symmetry in the 
structure of the ground states. 

4 Critical behavior and Local State Probability 

The model becomes critical as p — > 0. Exploiting the integrability of the system, it 
is possible to extract the critical exponents exactly. The inversion method [23] gives a 
method to calculate the critical exponent a. Unfortunately it turns out that the definition 
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Figure 5: Finite lattice with boundary heights (b m ,b m+ i) fixed in a ground state config- 
uration. 

of a is ambiguous for this model, since the free energy per site is regular as p goes to 
0, as was observed in |18j . Anyhow, in view of the form of the local state probability 
in equation (14.21) . there is only one natural definition for a consistent with the relation 
between critical exponents and conformal dimensions of the CFT : 



2-a = L (4.1) 

Using Baxter's transfer matrix method |23j, the Local State Probability (LSP) can be 
calculated exactly. The LSP P(a\b, c) is the probability that a given site of the lattice has 
height a, with fixed boundary conditions in the ground state (b, c). In order to calculate 
this quantity, we consider initially that the lattice is finite, with a shape as in Fig. [5], 

and fix the boundary heights to have the value they would assume in a particular ground 

{b if m is odd 
Finally, we take the limit m — > oo 
c if m is even 

when the lattice becomes infinitely large, all boundary sites being infinitely far from the 
center site. 



11 



Since this model is equivalent to a 4 fusion of the ABF model, the LSP has already 
been calculated in [18]. In appendix Owe rederive it in this particular case. The exact 
expression for the LSP is given by the following expression : 



7r(3 + 7]) (A + 1) \ . / ixdr \ . fKas\ 



siu i I sin I — I sru 



"J X{A,r;s}(P ) (4-2) 



where c = 6 + 77, ci = (6 + c)/2 — 1 and X{A,r ;s } are the branching coefficients of the 

coset SU fu(2) S k +l 2h ( cf appendix EJ) with < A < 4, 1 < r < Jfe + 1, l<s<k + 5. 

(2) 

The appearance of the branching coefficients of the Z3 (k) parafermionic theory in 
the local state probability is quite interesting. The branching coefficients encodes the 
conformal data, and they should a priori be relevant for the model only at criticality. 
But here these objects apply to the model away from criticality. However surprising, this 
is a standard result for this general class of RSOS models in regime III [21 J . 

A lot of informations about the conformal field theory can be extracted from these 
branching coefficients. In particular the central charge and all primary conformal dimen- 
sions are accessible through the limit q — > : 



The expression ( 14.21) is valid for any value of the nome p. The probability that the 
central site of the lattice has height a depends on the boundary conditions (b, c), even in 
the thermodynamic limit. This means that there is long range order in this model away 
from criticality (for < p < 1). The regime III of this model is ordered. 

Let us examine the expansion of the LSP around p = 0. Since the branching coeffi- 
cients behave as X{A,r- s }(t) ~ t~ c/24+Ar - s+5 (l + 0(t)), the LSP is dominated as p — > by 
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the branching coefficient of the identity operator (A = 0). At p = we find that the 
LSP becomes independant of the boundary condition {b, b + rf) : 

P(a\b 1 b + r ] ) = P° a = I sin[ T ) (4.3) 

Long range order disappears at p = : this is the signal of criticality. Moreover one 
can extract critical exponents by looking at the powers of p in the LSP. The presence of 
branching coefficients in the LSP (14. 2 j) forces critical exponents to be of the form : 



(3 = LA (4.4) 

where A is a generic conformal dimension of the parafermionic theory Zg (fc). a is 
taken as 2 — a = L in order to recover the standard relation between critical exponents 
and conformal dimensions : 



This value of a determines immediatly the conformal dimension of the energy operator 
through : 

a = ^ (4.6) 
2 - 2A £ v ; 

This leads to A e = 1 — and it corresponds to the following neutral descendant of 
the doublet T>i 3 : 



$e = V>L|^(1,3) = ^-1^1,3) ( 4 - 7 ) 

(2) 

It is interesting to note that the pertubation of the second parafermionic theory 1* 3 (k) 
by this field V'L2^ , (i,3) has been studied in the continuouum limit [251 12S]. In particular 
it is known to be integrable in the continuum case : an infinite set of integrals of motion 
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are preserved by this perturbation. The solvable model realizes this perturbation on the 
lattice. 

5 Discussions 

In this article solvable lattice models realizing at criticality the second parafermionic 
theory Z 3 ! (k) are studied. A version with positive Boltzmann weights is presented, for 
the purpose of eventual numerical applications. The critical behavior is described by the 
parafermionic zf\k) conformal field theory, as can be seen through the evaluation of 
the Local State Probability. 

Away from criticality these models describe the second parafermionic theory Z 3 (k) 
perturbed by the relevant field & e = ty^^Dx^ with conformal dimension A e = 1 — 2/(k + 
6). This perturbation by a neutral field should preserve the Z 3 symmetry, but it remains 
hidden in both the Boltzmann weights and the ground states. 

This is very different from the many solvable lattice models available for the first 
parafermions [3]. For instance the N state spin model [I], or the ABF model in regime II 
[22], in which the Ztv symmetry is explicit. The N state spin model has a Ztv invariant 
hamiltonian, while the ABF model in regime II realizes this symmetry at the level of 
ground states. 

(2) 

The model presented here for the Z 3 (k) parafermions is not new, however for the 
general case N > 5 no lattice model is yet available. I hope to discuss this point in a 
future publication. 

Acknowledgements: Very useful discussions with VI. S. Dotsenko are gratefully 
acknowledged. 
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A Explicit expression for the Boltzmann weights 



Explicitly the weights are given by the following formulas. Some Boltzmann weights 
(IA.10IA.14|A.18|) will be negative for u = —1/4, but it should be stressed that these 



signs are irrelevant and can be changed by a gauge transformation (cf section [2]). The 
weights are parametrized in terms of the elliptic theta function : 

7T7/ 

M = di( T ) (A.i) 

oo 

0i (u, p) = 2p 1 / 8 sin(w) Y[ (1 - 2p fe cos2u + p 2k ) (l - p k ) (A.2) 

k=l 

which enjoys the following properties : 

[L-u] = [u] (A.3) 
[-u] = -[«] (A.4) 

The following notations are used : 



H = ±1 (A.5) 

a M = a if jj = 1 

-aif/i = -l (A.6) 

A o( x) = ( Il + ^r a ' )' /2 (A7) 

For all the 14 weights given here \x can take the values ±1. All remaining weights can 
be recovered from reflexion and rotational symmetry. 



a a + 2fx\ _^ [1/2 + u][l + u] [3/2 + u] [2 + u] 

a + 2/x a + 4/i 

[1/2] [1] [3/2] [2] 

, a a + ^ \ = [1/2 + u] [1 + tj [3/2 + tj K + 3/2 + M ] 
a + fi a + 3/x 

[1/2] [1] [3/2] [o„ + 3/2] 
15 



w 



w 



w 



w 



w 



w 



a a + n 

a + 2/i a + 3/i 

a a + 2/i 

a + 2/i a + 2/x 



a a + /i 

a + /i a + 2/i 



a + 2/i a + 2/i 



a + /i a + 2/i 



a a + 2/i 

a + 2/i a + /i 



[1/2 +_u] [1 + m] [3/2 + u] [u] 
[1/2] [1] [3/2] [2] 



A 2 (a M + 3/2) (A. 10) 



[1/2 + u] [1 + u] [a M + 1 - u] [a M + 3/2 



[1/2] [1] k + 1] K + 3/2] 



[1/2 + u] [1 + u] f[l-u] [3/2 + u] 



(A.11) 



+ 



[1/2] [1] V [1] [3/2] 
M [1/2 + «] [a M -l/2][a M + 5/2] 



[2] [1/2] [a M + l/2][a M + 3/2] 
[1/2 + u] [1 + u] [1/2 -u] [u 



(A.12) 



[1/2] [1] [3/2] [2] 
A!(a M + 3/2)A 2 (a„ + l) 



Ai(a M + l/2) 



[1/2 + u] [1 + tj M [a M + 1 + tj 
[1/2] [1] [1/2] [a, + l] 

/[1/2][1] [a, + 5/2][a M -l]\ 1/2 



V[3/2][2] [a^ + 3/21^ + 1], 

[1/2 + tj [a„ + l-tt] [a^ + l/2-tt] 
[1/2] [a, + 1] [a,, + 1/2] 

[ty + 3/2 - u] 
K + 3/2] 



(A.13) 



(A.14) 



(A.15) 



W 



W 



a — fi 
a + ji 



a 

a + /i 



[1/2 + tj [tj [a„ + tj [ty + 1/2 - tj 
[1/2] [2][a„-l/2] [a,. + 1/2] 

/[3/2][2]y /2 / [a, + 3/2] [a,- 3/2] \ 1/2 
V[1][1/2]J V [«m + 1]M J 

[1/2 + u] [cy + 1/2 + u] ( [1/2 -«][! + u 



[1/2] [1] 



+ 



[1/2] [a„ + 1/2] 

M[l/2 + ^] [3/2] [ a ^-l][g^ + 3/2] \ 
[1/2]* [2] [a,+ l][a,-l/2}J 



(A.16) 



(A.17) 
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a + fi a + /i 
[1/2 -u] [1+u] 
[3/2] [1/2] 




+ 



[u] [1/2 + u] ( [a„ - 1] [a M ] [a„ + 1] [a M + 2] \ 1/2 



l/2][ aAt + l/2]2[o M + 3/2]J 



[1/2] 1 



N [1/2 + u] [3/2] [a M -l/2][a M + 3/2] 



[1] [1/2] [2] 



M [a^ + 1] 



[a M - w] [a M + 1/2 - w] [a^ + 1 - it] 
[a M ] [a„ + 1/2] [a„ + 1] 

[ty + 3/2 - u] 

K + 3 / 2 ] 



(A.18) 



(A.19) 



W 



a 



a + fi\ _ [a M - n] [^ + 1/2 - u] f [1 + n] [1/2 - n] / [3/2] 2 [l] 



a + /i a 



[a, + 1/2] 1 [1] [1/2] l[l/2H2] 
[2]^ + 1/2] 2 \[u][l/2 + u 



a a 



a / i a M 



[l][a M -l/2][a M + 3/2]; [1] [1/2] 

[a M - it] [a M + 1/2 - tt] [a^ + 1 + m] [a M + 1/2 + u] 
[a M + 1/2] [a M + 1] [a„ + 1/2] 



(A.20) 



+ 



M [l/2+jj [a^-3/2][a M + l] 
[1] [1/2] [a,-l/2] 2 
) [a^-u]{a^ + 1/2 + u] [l] 3 
1 [a M + l/2] 2 [l/2] 2 [2] 
[l + u][l/2-u] [gg + 1/2] [a M - 1] [a M - 2] 
[2] [3/2] [a, + l] 2 M 

(2), 



(A.21) 



B The second parafermionic theory Zg (k) 

The details of the second parafermionic theories "L^\k) can be found in |3]. These 
CFT possess an infinite symmetry generated by parafermionic currents. The Virasoro 
generators form a subalgebra of this enhanced algebra. The second parafermionic theories 
are realized by the following coset construction : 



bU {Z) k+4: 



and possess a central charge : 
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-Ml - „ , " , oN ) (B.2) 



(A; + 2)(fc + 6) 

The chiral algebra is made of 2 parafermionic currents *ff, of dimension 4/3, and 
Z3 charge ±1. They obey the following operator product expansion : 

x * -> tft (B.3) 
* x tft _^ J (B.4) 

The parafermionic algebra primaries of the second Z 3 conformal theory are labeled 
by 2 integers (r, s), with conformal dimension : 



_ ( r(fc + 6 )- g (fc + 2)) 2 -16 m(4-m) 

r ' s ~ 16(k + 2)(k + 6) + 48 1 ' 

m = s — r mod 4 (B.6) 

1 < r < jfe + 1 

1 < s < k + 5 



The value of m labels different sectors of the theory. Singlets (m = 0) are neutral 
fields, while doublets (m = 2) have a Z 3 charge ±1. Finally m — 1,3 correspond to 
disorder fields. 

In each sector (r, s) , descendants will have a conformal dimension A r s + 5 + N de- 
pending on their Z3 charge (cf Fig. EJ 



5 = or 1/3 in a singlet module 
5 = or 2/3 in a doublet module 
5 = or 1/2 in a disorder module 
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q = 



q = +l 




-1/3 1 



1/3 



*_1 /3 <S 



O 



o 



Figure 6: Structure of a singlet module. Arrows depict the first (charged) descendants 
of the primary singlet S, obtained by action of the parfermionic mode operators. 



q = -l 



-Q- 



q = 



q = +l 
V 




2/3 



Figure 7: Structure of a doublet module. Arrows depict the first (neutral) descendants 
of the primary doublet T>,T>\ obtained by action of the parfermionic mode operators. 
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In particular the conformal dimension of the neutral field e = V'La-^Ci^) = V'-^^tl 



(1,3) 



k + 4 l-T (B.7) 



k + 6 L 

SU(2) k xSU(2) 
SU(2) k+4 

the LSP. They are expressed in terms of level 4 string functions of SU (2) c; 



The branching coefficients of the corresponding coset — su(2) k+ 4 appear naturally 



X{ A,r;s} = Tr {A , r;s} (q L °- c ^) = ciF m {q) (B.8) 



m=0 



F m (q) = Yl 5 ™,miq Pr - 2 ^ k+6) - a - (W^— 2 "^ 6 ^ (B.9) 

n=— oo 

and 

_ ((fc + 6)r-(fc + 2)s) 2 
^ " 16(lfc + 2)(fc + 6) (B - 10) 



m! = r — 2n(k + 6) — s mod 8 
m" = — r — 2n(A; + 6) — s mod 8 

In particular as q — > they behave as : 

X{A,r;s}(q) ~ ^ s+<JA , r>s - C /24 (l + ^ (Rn) 

where <5A,r,s depends on the subsector of the primary field $ riS : 

• if r — s = mod 4 : then m! = 0, 4 mod 8. Using = c^I^ we can always take 
m' = 

Then we get : 

q 2/24 4 ~ 1 
20 



g 2/24 c 2 „ g l/3 
2/24 4 1 

reproducing the structure of the module of a singlet : neutral descendants of a 
singlet have a conformal dimension in A r>s + N, while charged descendants have a 
dimension in A rjS + 1/3 + N 

• if r — s = 2 mod 4 : then m' = 2, 6 mod 8. Using = c?_ m we can always take 
m' = 2 

In this case we get : 

q 2 ' 24 cl ~ g 1 / 12 
g 2 / 24 c° ~ gVl2+2/3 

which is precisely the structure of the module of a doublet : charged descendants 
correspond to S = 0, and neutral descendants to 5 = 2/3. 

• if r — s = ±1 mod 4 : then m! = 1, 3, 5, 7 mod 8. We can always take m! — 1 

g 2 / 24 cj ~ g 1 / 16 

g 2/24 C? ~ gVie+l/2 

describing the disorder sector. 

C Local state probability 

Since this model is strictly equivalent to a 4 fusion of the ABF model, the LSP has already 
been calculated in [IB]. In this appendix this result is rederived in a slightly different 
manner. In particular a simpler form is given for the non restricted one dimensional sum, 
more in the spirit of [2~T] . 
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Using the corner transfer matrix method [23J, one obtains for the LSP : 

P(a\b,c)= lim P m (a,b m ,b m+1 ) (C.l) 



with 



P m (a,b,c) = Sj-UaXmta, b, c;x 2X 

u a = E{x a ,x L ) 

S m = y^M a X m (q,6, c) 

X m (a,b,c;q) = 



(C.2) 



a,2,...,a rs 



where the sum is taken on all admissible paths (02, ...a m ) from a% — a to a m+ i = b m +\. 
The problem is now to evaluate the one- dimensional configuration sum X m , and then the 
normalization S m . 

The restricted ID sum X m (a\b, c; q) is uniquely determined by the recurrence : 



X m (a,6,c;g) = ^ X m-i(a, a m , b)q mH 



(a m ,6,c) 



'-tin 



X (a,b,c;q) = S a)b 



(C.3) 
(C.4) 



where the sum Yl" is taken over a m such that the pair (a m , b) is admissible, and 
if (a, b, c) — H(b — a,c — a) has the following form : 





-2 


-1 





+ 1 


+2 


-2 


2 


3/2 


1 


1/2 





-1 


3/2 


1 


1/2 





1/2 





1 


1/2 





1/2 


1 


+ 1 


1/2 





1/2 


1 


3/2 


+2 





1/2 


1 


3/2 


2 



(C.5) 



Let's consider initially the non restricted one- dimensional configuration sum : 

f m (a,b,c;q) = Y,'frn-x{a,a m ,b)q mH ^ b ^ (C.6) 
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f (a, 6, c; q) = 5 a , b 



(C.7) 



where the sum is now taken over a m such that the pair \a m — b\ = ±2, ±1,0 (i.e. 
we forget the bounds over the heights). This quantity is invariant under translation, so 
that it depends only on (7, 77) = (b — a,c — b) 



The solution can be expressed in terms of q-multinomial coefficients 



m 
k\...k r . 



h =m, ki> 



with 



i=l 



The unrestricted ID sum is then of the form : 



m 



(C.8) 



(C.9) 



k 

where k = (k_ 2 , fc , fci, k 2 ), the sum ^£ is taken over k such that : 



m 
ki 



(CIO) 



fc,- > 



i=-2 
2 

^2 iki 



= m 



7 



i=-2 



and <2(/c) is the following quantity : 



Q(k) = Hh - 1) - 2A; 2 A;_ 2 + (k 2 - k- 2 )(h - fc_i) + k {h + fc_x) (C.ll) 



This can be simplified using the constraints on fc: 
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Q(k) = 7 2 /4 -m+(k 0+ *I±*=l) 2 + ^±^i (C.12) 

The restricted ID sum is then the usual sum over all reflexions along the boundaries. 
We get : 

X m (a,b,c;q) = Y m (a,b, c; q) - q~ a Y m (-a,b, c; q) (C.13) 

oo 

Y m (a,b,c;q) = £ q~ 2 ^ L+ka ^ +kL f m (b - a - 2kL, b - c; q) (C.14) 

fc=— oo 

It is straightforward to check that this formula for X m obey : 

X m (a,0,2;q) = 
X m {a,L,L-2;q) = 

The other cases are more involved. For rj e {±2, ±1,0} the following equations are 
satisfied by X m : 

q mH ^X m (a,l/2,3/2;q) + q mH ^X m (a,-l/2,3/2;q) = 

g m ^X m (a,l,l;g) + g m ^X m (a,0,l;g) + g m ^X rn (a,-l,l;g) = 
q mH(-i, v)Xm ^ g/^ 1/2 . g) + ^(o,,) Im((() 1/2? 1/2 . g) + 

g m ^X m (o,-l/2 ! l/2;g) + g m ^X m (a,-3/2,l/2;g) = 

All this relations ensure that the restricted recurrences are satisfied by X m . 

The limit m — > oo of of f m (l, f]\ q) has the following form can be taken by expressing 
hi-, k-2 in terms of the other variables ki,ko,k-i : 



Ak 2 
4A;_ 2 



= 7 + 2m - (3fci + 2A; + fc_i) 
= -7 + 2m - (fci + 2A; + 3fc_i) 
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The sum is now on [ki, k , k-x) such that k 2 and k_ 2 are integers : 2m + 2k = 
k\ — k-i + 7 mod 4 

Moreover the rj dependent term J7 kiH(i,T]) takes the form : 



m + 7/2 if 1] = 2 
3fci + 2k + fc_! 

+ 7/4 if 7? = 1 

fci + 2/c + 

m it rj = 

2 ' 

fci + 2/co + 3Ar_i ., 

- 7/4 g lf 77 = - 1 

m — 7/2 if 77 = —2 



hm f m (l,v) = q^' A -±- T 



1 



,, h+ ,./ 4 _ t(so „ )( , ) 



where Q v {k) is given by : 



Ko H ^ I + 



3fci+2fc +fc-i 



k 1 +2k +k. 



fci+2fc +3fc_ 



2 





if 77 = 2 

if 77 = 1 

if 77 = 

if 1] = —1 

if r] = -2 



(C.15) 



(C.16) 



* stands for the restriction 2/co + £4 — k-\ = Sq mod 4, and So = 2m — 7 mod 4. 

Up to a multiplicative factor, the functions $( SOjJ) )(g) are precisely the branching 
functions of the c = 1 Z 4 parafermions [2T] : 



A, ,(V) - fl l/24-(2+r ? )(2-r,)/48,2+^ 



(C.17) 
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The branching functions b l m are defined for / mod 4, and m mod 8, I — m even, and 
enjoy the following symmetries : 



bl 



b l - b l 



li _ iA-l 

°m — °4-m 



Finally we get for the limit of f m : 



lim f m (l,v) = g^-^)/ 4 -^gV 2 4-( 2+ ,)(2-,)/48 C;? _ 2 

?Kry /24 - (4 -^ )/48 C V, 



,7(7+r;)/4 



= Q 

= ? 7(7+J?)/4 / /48 c^ l 2 _ r) (C.18) 

where c l m = ri(q)~ l b l m is a level 4 string function of SU (2). We now impose the ground 
state (6, 77) at the boundary : 



b if m is odd 

(C.19) 

c = b + i] if m is even 

and compute the thermodynamic limit m — > 00 of f m (b m — a, b m+ i — b m ; q) (in the 
case m odd for instance, but the result does not depend on the way the limit is taken ) : 

Let's say that in the case m odd we have a certain s = 2 — (b — a) mod 4. Then for 
m even we will have s' = — {b + r\ — a) = sq — (2 + rj). 

lim f m (b m - a, b rn+1 - b m , q) = q^ b+ ^ A q^ ^c^_ 2 _ v if m is odd 

TO— >00 

lim /TO(&TO-a,Wi-&TO;?)=? (6+ ^ a)(6 - a)/4 / /48 47-2-, +4 if m is even 
But the symmetries of the branching coefficients (b l m = b^l A ) ensure that c 2 ~'L 2 _^_ | _4 = 

2+r? 
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lim f m (b m - a, b m+1 - b m , q) = q^ a+ ^ V /M c+ V (C.20) 

m^oo 

where m! = 2a - 2dmod8, d = (6 + c)/2 - 1 = 6 + 77/2 - 1 

Recalling that L = fc/2 + 3, we recognize in linim^oo X m (a, b m , b m +i) a SU< su(2) S k +^ 4 ' 
branching function (cf appendix [B| : 

hm X m (a,6 m ,6 m+ i) = g v «^- 2 > X{2+^,2d ; 2a} 

m^oo 

= ? (d+1 - a)2/4 - 2/24 - M ^^X {2+ , )M ; 2 a } (C21) 



As usual in regime III the principal specialization of the character identity arising 
from the coset construction comes very handy to normalize the LSP. 



(fe+4) 

P(a\b,v) = X{K2 ^ 2a - 1 (C22) 

X\ X2d-1 



^+1-^/4-^/24- ^Z^ 2 (a23) 

U ' r "" rL) X^W*) (C.24) 



E(x d } x L - 2 )E(x( 3 +v)/\ x z) 
where we introduced : 

00 

E(z,x) = Y[ (l - zx n - 1 ) (l - «~ 1 x ft ) (1 - x n ) 



n=l 
00 



£ 

fc=— 00 

00 

= II(i- 



n=l 

In order to study the critical behavior of the model, described by the LSP around 
p = 0, it is useful to go back to the variable p rather than x. Then again, having identified 
in the LSP a ratio of characters is of great help since modular transformations are readily 
available. Finally we get the following expression : 
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P(a\b,c) — v 



with : 



Vl\K L -2>P )°l\ n 6 > K- r -s=Q mod 2 

. /7r(3 + 77) (A + 1)\ . / ndr \ . /nas\ . r . nr . s 

8111 { 6 ) Sm \L^2j 8m [~TJ ^ a ^}(p ) ( C - 25 ) 



i] = b — c 



and the sum in the r.h.s. is over : 



< A < 4 

1 <r < 2L-5 = k + l 
1 <s< 2L-l = k + 5 
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